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Abstract. We continue an earlier study of the shape parameter c contained in the famous mul- 
tiquadrics (-l)^(c 2 + \\x\\ 2 f , /3 > 0, and the inverse multiquadrics (c 2 + ||a;|| 2 )' 3 , f3 < 0. In [5] 
the space of interpolated functions consists of bandlimitcd functions. Now we are going to treat a 
more general function space which roughly speaking is the same as the native space of gaussians. A 
totally different set of criteria for the optimal choice of c will be provided. 
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1 Introduction 

As before, we are going to adopt a seemingly more complicated definition 

h(x) := T(-^)(c 2 + \x\ 2 )?, (3 e R\2N> , c > (1) 

, where \x\ is the Euclidean norm of x in R n , T is the classical gamma function, and c, /3 are 
constants. This definition will relieve our pain of manipulating its Fourier transform and developing 
useful criteria. 

Recall that h(x) is conditionally positive definite(c.p.d.) of order m = max{[^], 0} where [~|] 
denotes the smallest integer greater than or equal to ^. This will be used in the text. 

For the reader's convenience we review some basic features of the development in [5]. For any 
interpolated function /, the interpolating function will be of the form 

N 

s(x) := Cj h(x - Xj) + p(x) (2) 
»=i 

where p(x) € P m -i, the space of polynomials of degree less than or equal to m — 1 in R n , X = 
{xi, • • • , xn} is the set of centers(interpolation points). For m = 0, P m -i '■= {0}. We require that 
s(-) interpolate /(•) at data points (xi , f (xi)) , • • ■ , (x n , f (x n)) ■ This leads to a linear system of the 



1 



form 



N Q 

^2 (kh(xj - Xi) + ^2b i p i (x :j ) = f(xj) , j = 1,- • • , N 



(3) 



N 



^CiPjixi) = , j = 1, • • • ,Q 



to be solved, where {pi, ■ ■ ■ ,Pq} is a basis of P m -\. 

The sovability of the linear system is guaranteed by the c.p.d. property of h. However if c is very 
large, h will be numerically constant, making the linear system numerically unsolvable. Moreover, as 
pointed out by Madych in [8], if c is very large, the coefficient matrix of the linear system will have a 
very large condition number, making the interpolating function s unreliable when f{x%), • • • , /(xjv) 
are not accurately evaluated. 

Each function of the form (1) induces a function space called native space denoted by 
Ch.m(R- n ), abbreviated as Ch, m , where m denotes its order of conditional positive definiteness. For 
its definition and characterization we refer the reader to [2] , [3] , [5] , [7] and [5]. This space is closely 
related to our space of interpolated functions. 

As in [5j , we need the following basic definitions for our development of the criteria. 

Definition 1.1 For n = 1,2,3, the sequence of integers "/ n is defined by 71 —2 and 7„ = 
2n(l + 7„_i) ifn>l. 

Definition 1.2 Let n and (3 be as in (1). The numbers p and Ao are defined as follows. 

(a) Suppose /3 < n — 3. Let s — |~ "~^~ 3 ] . Then 
(i) iff} < 0, p = 2J- and A = Q±^t± ; 

(ii) iff»0, p = l + ^ and Ao = 
where m = \^~\ . 

(b ) Suppose n — 3 < < n — 1. Then p — 1 and Ao = 1 . 

(c) Suppose $ > n - 1. Let s = - \ n ~%~ 3 ] . Then 

p = 1 and Ao = 7- -77- -r 7- -r where m = {^-] . 

H (2m + 2)(2m+l)---(2m-s + 3) '2 1 



The following theorem is cited directly from [5]. 

Theorem 1.3 Let h be defined as in (1) and m = max{0, [§]}■ Then given any positive number 
bo, there are positive constants Sq and A, < A < f , which depend completely on bo and h for which 
the following is true: For any cube E in R n of side length bo, if f € Ch, m and s is the map defined 
as in (2) which interpolates f on a finite subset X of E, then 

\f(x) - s(x)\ < v^Tc* v^(A)* ll/IU (4) 
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holds for all < 6 < Sq and all x in E provided that 8 — d(E,X) :— sup^g^ inf^gjf \y — x\. Here, 
a n denotes the volume of the unit ball in R n , and c, Aq were defined in (1) and Definitionl.2 
respectively. Moreover Sq — 6( g 7 \ m +i) > an d ^ = (|) eo ~ ln where 

2 1 ,9 



C = max|2p'^e 2 ^", J-| , p 



c 



The integer 7„ was defined in Definition 1.1, and \\f\\h is the h-norm of f in C/i )fn . The constant p 
was defined in Definitionl.2. 

Remark: Obviously the domain E in Theorem!. 3 can be extended to a more general set fi C R n 
which can be expressed as the union of rotations and translations of a fixed cube of side bo. 

In this paper the space of interpolated functions is defined as follows. 
Definition 1.4 For any positive number a, 



E a := {/ e L 2 (R n ) : J |/(£)| 2 e^ < oo} 



where f denotes the Fourier transform of f . 

Remark: It's easily seen that E a is just the well-known native space of gaussian. For each / in E a , 
we define its norm by 

ll/IU, ■■={! iKOfe^dt 



2 Fundamental Theory 

It's easily seen from Theoreml.3 that the error bound (4) is greatly influenced by the choice of c. 
This is indeed the starting point of our theory. However, in order to develop useful criteria for the 
choice of c, some technical manipulation and theoretical analysis are necessary. 

Lemma 2.1 Let a > and f3 < 0. // \n + f3\ > 1 and n + [3 + 1 > 0, then E a C C h , m {R n ) and for 
any f e E a , 



1/2 



kk£ I , . , n+p+i c |f|_JSl; 



\\f\\h<2- n - — Tr- n —*c^{ sup \^\^e c ^-—} \\f\\ E „ 

[\t\eR« J 

where \\f\\h is the h-norm of f in the native space Ch, m {R n )- 
Proof. Let / G E a . By [7] and [2], 

1/2 



if ii i 1 f^'V 

l/lu - WFl~W^\ 



(27T) 



1/2 

l/(C)| 2 



2^ (f)" 2 " 7 x:^(c|ei) 



d^ ) {Theorem?,. .15 of [9]) 



1/2 



< 



(27T) 



l/(OI 2 V^|e^l 



-d£ > (Corollaryh.Yl of [9]) 



l-n-<3 

C -1 



2 n+ 



< 



n 1/2 

l C l^ e ci?n 1/2 

, , ) SUp [> \\f\\ Ea 



1-Tl-ff 

C 4 



< oo. 

The lemma thus follows immediately by Corollary 3. 3 of [?]■ 

In the preceding proof we didn't find the supremum of a function. Let's try it. Suppose 

G(0 ^r^e^f £>0 



. Then 



G'(0 = 



-r 



r 



c--e 

(7 



iff 



iff 



. So, 



z fT 



= ccr + y/cV 2 + 4a(n + [3 + 1) =< 



< 2 



,l 4 7T n 4 C 4 



(O 



2 CT 



11/11 



. We sum it up in the following theorem. 

Theorem 2.2 Under the conditions of Lemma2.1, any f £ E a satisfies 



i+j3 l l-n-p 

h < 2- n - — TT- n -- C —^- 



(O 



„t* err 
2 e c « 



11/11 



where 



e 



+ ^/c 2 o- 2 +4a{n + f3 + I) 
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Corollary 2.3 Let a > and /3 < 0. If \n + f3\ > 1 and ?i + /3 + 1 > 0, f/ien ^ in Theoreml.3 has 
the form 



\f(x) - s(x)\ < 2-4"7r-4"7^VA c— 
whenever } € hj a , where t, ■= - r . 



...... '» + .«+! „t* 

(f) 2 e c « 



(A)*||/|| 



Remark: Note that Corollary2.3 covers the most useful cases (3 = — 1 and n > 2. However the case 
P = — 1 and n = 1 is excluded. For /3 = — 1 and n=lwe need a different approach. 

Lemma 2.4 Lef cr > 0, (3 — —1 and n = 1. Then E a C Ch, m (R n ) and for any f E E a , 

t \ - 



\\f\\H< 



2«+3 7T 1 ^2 

Proof. Let / e E a . By and 0, 



' 2V§ sup v^e ^!" 1 ^ 



iei>i 



11/11 



II /Ik = 



1 



l/(0P 



d£ 



l/W 



(2tt)» U V2/C (c|e|) 



d£ > (17ieorem8.15 0/ [9]) 



l/(OI S 



\m\ 2 



Here. 



I/(Q| 2 

i £ i<i Mc|£l) 



dcS, 



(2^)«2i \J mi KoWt\) J m> ± KoMt\) 

l/(0l 2 



|f|<i -ln^ 



d£, (fco( z ) ~ — {(^g) r J'^°( z ) as z — > where r ~ 0.577) 



oy p.255, p.374 and p.379 0/ [T] 

lei 2 1 



l«l<7 



< sup 



l/(0|V 



UH 1 , c|£| 
e " l/n-j 1 



l/(e)i 2 e^^ 



< 



< 



e " m-^-M 



|L ^ere0< |r| < 



|Zn2| 



ll/l 



. Also, since = y/n, 



f Wm^ < 2^/ 1/(01 
J\e\>i ^(cici) /|f|>i 



=ciei 



2 nil Vc|e|e c 

e - jjg— 

e <* 



< 2V3sup ^|^-||/|||, 
l«l>! e~ 
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. Our lemma thus follows immediately. jj 
Theorem 2.5 Let a > 0, f3 = — 1 and n = 1. -For any f e i^, 

11/11, < |^ + 2V3M(c) J 2 ||/||^ 

w/iere M(c) := e 1- ^ i/ c < ^ and M(c) := g ^+>W+4 g^ ^ c > _2_^ w/lere g (£) := 

Proof. The maximum of on [^,oo) obviously exists. In order to find its exact value, we first 
find g'(£). Note that 



g'(0 



e " 



Hi 
e a 



2^/1 V S CT V S 



Z 2 

e " 



iff 

iff 
iff 
iff 
iff 



Vc + 2cVc£- -\/^ 2 =0 



1 + 2c£ - -<f = 



CO" + \Jc?a 2 + 4a 



Let ^ , = cg +vc^+4 ff - Then £* < I iff c < Also, lim^ + .9(0 = lim^oo - 0. This 

gives that 



sup fif(0 



> 



6 



. The theorem then follows from Lemma2.4. ft 

Corollary 2.6 Let a > 0, j3 — —1 and n — 1. For any / G formula (4-) in Theoreml.3 can be 
expressed as 

i 

\f(x)- S (x)\<2^n'^^ l ^(X)^l=^ + 2V3M(c)y \\f\\ E „ 
Now we begin the study of the case j3 > 0. 

Lemma 2.7 Let a > 0, /3 > and n > 1. TTien C Ch, m (R n ) and for any f 6 



l-g-r. 



ll/IU < rfoc^ ^ sup = ||/|| Ba 



where d Q is a constant depending onn,/3 only. 
Proof. By definition, 



|a|— m 



ii/ii" = { E ^\\( Da f) A \\U P )\ m°>ndm 



E 5/ \( Da n A (o\ 2 d P 

y\—m 

Y 4 / i(^ Q /) A (£)i 2 — = — <«e J- 

y 4 / |i m r/(oi 2 1 -^—dt 

M)* { v i [ e a \Ko\ 2 d t 

n\ I \C\2mh(C\ 



(27r)"V2 1+ # [ ■/ c 2 £«££(c|f|) 

where C(m,n) denotes the number of terms in 



(2^)«V2 I +l r c £ ? 1 



7r e~ 
2 



Vm!C(m,n) ^ ( A fl\ f [ \f(^\^ 
(2^)«V / 2 1 +4 
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< d c 



sup 



lei 2 
e o 



where do :- 



y / m!C(? 



Since \\f\\ h < oo, / e C h)T 



Theorem 2.8 Under the conditions of Lemma2.7, 

||/IU<doc^§- > 



«*) 2 
e =• 



* _ cg+^/e 2 <T 2 +4g(l+ff+n) 
4 



w/iere £ 
Proof. Let g(x) : 



l + g + n 

z ^7 eCX , x > 0. Then 



2 [ „2 



g'(x) = e - < e 



iff 
iff 

iff 
iff 

iff 



1 + [3 + n 



x 2 e + x 2 e c 



- a; 2 e e 



= e " e 

ic 2 

= e c "— 
= 



1 + (3 + n + l + g + n l + g + n 2 

x 2 + cx 2 — a; 2 — x 

2 (j 



~x 2 



1 + /3 + n _i i i 2 
a; 2 + cx 2 — x 2 — x 

2 (7 



l + /3 + n 12 

• — ^ + c\/x iV 1 = 

2 yjx a 



4 „,2 



1 + 13 + n + 2cx - -x 
V5 = U 

Ax 2 - 2cax - cr(l + p + n) = 



£ = 



2ca + y/Ac 2 cr 2 + 16cr(l + (3 + n) 



co- + ^c 2 a 2 + 4(7(1 + /3 + ■ 



. Then the theorem follows immediately from the preceding lemma. 



Corollary 2.9 Let a > 0, (3 > and n > 1. for any / e £7 CT; (^J m Theoreml.3 can be expressed 
as 

\f(x) s(x)\ < 2^^^^X~ {X)U c 1±i ^ | (g!2lll^ | ||/||^ 
w/iere do *s defined as in Lemma2. 7, and £* is defined as in Theorem2. 8. 



8 



3 How to choose c?— a more practical approach 



Our criteria for the optimal value of c are based on the exponential-type error bounds introduced 
in the preceding section. Intuitively, the core of those error bounds is (A) s which converges to zero 
as the fill-distance 5 tends to zero. Surprisingly, this is wrong. In practice (A) 7 can be essentially 
ignored, as explained in [5]. What's influential is the other part determined by c. In this section we 
will ignore (A) « totally and develop concrete criteria for the optimal choice of c. 



Casel. /3 = — 1 and n > 2 Let / G E a and E be the cube in Theoreml.3 with side length bo- 
Let h be the map defined in (I) with (3 = —1 and n > 2. For any fixed 5 in Theoreml.3, if 
< S < 47 (m+i) ' tbe optimal choice of c is c* G [12/9^ne 2 ™ 7 "7„(m + 1)5, oo) such that i?(c*) is the 
minimum value of H(c) in [12p^/ne 7n 7 n (m + 1)<5, oo) where 



H(c) := c~% 



4n 

c+\lc z H 



(5) 



The constants 7„ and p were defined in Definitionl.l and Definitionl.2 respectively. 



Reason: Note that in Theoreml.3 there is a restriction 5 < So- That's why we put the restriction 
S < 4 7n (m+i) - The number 4 7 „(m+i) is J ust ^ obtained b Y letting C = ^ when c = Sfro/V^e 2 ™ 7 " ■ 
Increasing c(i.e. c > S&oPv 7 ™ 62 ™ 7 ™ ) does not change So, but decreasing c(i.e. c < ZboP\fne lnin ) makes 
So smaller. After S, where < S < 47 feqrn > * s fi xe d, c cannot be less than 12y0i/ne 2 ™ 7 "7„(m + 1)5 
due to the restriction S < So- 

When (3 — — 1 and n > 2, the crucial part of the error bound in Corollary 2. 3 influenced by c is 



c 4 



C 4 4 4 (7 * 





2 c 2 CT + cv / c 2 cr 2 +4crTl 




e * 



_1_ (ca- + y / c 2 o- 2 +4crTi) 2 
' 16 



c 2 + ^ 



■ [c 2 +cV^T^] 





g3 1 T [c 2 (T 2 +2c(Tv / c 2 (T 2 +4cm+c 2 cr 2 +4crn] 




4 cr 4 










- [c 2 ct 2 +c<7\J c 2 cr 2 +4crn+2(7nl 


4~?cT 


"f erf 


c + V c2 + v 


4 e f ^ 2 W- 2 +^J 



3 ^[c 2 + C ^+¥+ir] 



. We thus define 



•= r. 4 



4n 



The minimum value of H{c) on [12py / ne 2 ™ 7 ™7 n (m-|-l)(5, oo) is then what we want. 



Note that in the preceding criterion there is an unpleasant restriction c G [12p-\/ne 7 ™7 n (m + 
1)5, oo). Our fundamental theory (Theoreml.3) is not strong enough to replace it by c G (0, oo). This 



9 



is a question and deserves future research. However, decreasing S will decrease 12p^/ne 2nin ') n {m + 
1)6. The smaller 5 is, the more meaningful our criterion will be. 

Also, it's interesting and useful to know whether the optimal choice c* is unique. The answer 
is the affirmative. Let's show it as follows. 

By simple calculation, 



H(c) 





/ An 


1 + 









I e s 



. Therefore, 



H'(c) = 



n 



' + 1 


h 


An 

C 2 (J 


+ 1 


A + 


An 


c 2 a 



-An 



16 



2c- 



4n c 2c 



c z + — + - 



= e 







/ 4n 


e s 


1 + 











-An 



<?a x \ 



4n 





I An 


(7 




+ 


1 + v 1 + — 




2c+ \l 




V C Z <7 


16 





An 



c 2 + ^i 



iff 



iff 



iff 



-An 



1 + 



4n 





A 4n 


(7 




1 + 






2c+ \l 




V tr<7 


16 





CTC 3 W1 



4n 



+ 



An 



An 



a 
16 



2c- 



4n 



16 



/ 9 471 


/ 9 471 








c+ v c + - 




2c + y 











4n 

'c^ + — + 



n 



(6) 



In the left side of the last formula there is 



F(c) := 



c 2 



. Let 



Vc 2 + fe 



where fc := ^. Then 



F'(c) = 



c 2 + k 



c 2 + 2k~ 



> 



10 



for all c > 0. Therefore the left side of (6) is an increasing function of c. It follows that the optimal 
choice c* is the unique value of c satisfying H'(c) = 0. We conclude that the optimal value c* is 
unique. 

Remark: By the definition of H(c), it's obvious that H(c) — > oo both as c — > + and c — » oo. 
Consequently the error bound will be very poor when c is too large or too small. 

Numerical Result: The optimal value c* of c is theoretically where the minimum value of H (c) 
occurs. The value of c* can be obtained numerically by Mathematica or Matlab in a very easy way. 
One needs only find the number minimizing H (c) or solve the equation (6). Both are straightforward 
and involve very simple commands. 

Note that in the error bound (4) of Theoreml.3 there is a main function determined by c. Let's 
call it MN function, and its graph MN curve , in honor of Professor W.R. Madych and Professor 
S.A. Nelson. In our current case the MN function is just H(c) of (5). Sometimes the MN function 
will be multiplied by a constant to make its graph look better. Then the function and its graph will 
be called modified MN function and modified MN curve respectively. 

In this paper all MN curves start from c = 12/9 V / ne 2 ™ 7 ™7„(TO 4- 1)6, unless otherwise stated. 
Let's call this number the original left endpoint. In this paper we sometimes start the MN curve 
at some point to the right of the original left endpoint to make the graph look better. These are 
Fig.22,23,24,25,26,31,32,33,37,38,39,40,47,and 48. However in any case the lowest point of the MN 
curve always corresponds to the optimal choice of c in the entire interval [12py/ne 2n ' r "j n (m + l)d, oo). 
Nothing is lost. Now let's see some examples. 



Graph of the MN curve with 5=10 



,-22 



MN(c) 



4x 10 7 



2x 10 7 



6x 10 7 



8x 10 7 




c 



10.5 



1 1.0 



1 1.5 



12.0 



Figure 1: Here n — 2, /3 



1,(7 = l,6o 



1 and 5 = 10" 
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Graph of the MN curve with 6=10 23 

MN(c) 




1.1 1.2 1.3 1.4 1.5 

Figure 2: Here n = 2,/3 = -1, a = 1, b = 1 and £ = 10 
Graph of the MN curve with <5=10 -24 

MN(c) 




0.5 1.0 1.5 2.0 

Figure 3: Here n = 2, j3 = — 1, a = 1, bo = 1 and <5 = 10 
Graph of the MN curve with (5=10 -25 

MN(c) 




Figure 4: Here n = 2, /3 = -1, a = 1, 6 = 1 and <5 = 10 
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Graph of the MN curve with 6=10 



MN(c) 
40 - 

30 

20 

10 



Figure 5: Here n = 2, /3 = — 1, a = 1, &o = 1 and S = 10 



The second case can now be introduced. 



Case2. 



(3 = — 1 and n = 1 Let f <E E a and be the cube in Theoreml.3 with side length bo- Let /i 
be the map defined in (1) with /3 = — 1 and n = 1. For any fixed 6 in Theoreml.3, if < S < 47 r^ffj > 
the optimal choice of c in the interval [12p v / ne 2 ™ 7 "7„(m + l)5, oo) is c* 6 [12 ( o- v /ne 2 ™ 7 ™7 n (m + l)5, oo) 
such that H(c*) is the minimum value of H (c) on the interval [12py / ne 2n7 "7„(m + 1)(5, oo) where 



H(c) 



l 
1 

n/5 



I ^ + 2V3e 1 -^ 



7^2 + 2 v / 3G v 



ccr+Vc 2 a 2 +4:a 



if c< -|= 
if c> -2= 



, where 



G(0 



£1 

e " 



. The constants 7„ and p were defined in Definitionl.l and Definitionl.2 respectively. 

Reason: This is just a simple result of Corollary2.6. The reason we put the restriction c 6 
[12p v /ne 2n7 "7 n (m + 1)6, oo) is the same as Casel. ft 

Remark: Obviously H(c) — > oo both as c — > + and c — > oo. Therefore one should avoid let- 
ting c — > + or c — > oo. 

Numerical Result: Numerical experience shows that most time H(c) is increasing on 
and the minimum value of H(c) happens around c = -4=. In fact, on the interval (0, ) 

1 . rz , i_ /2 - cV 



3a 



H'(c) 



\ln2\c 2 



■2V3e 1_ 
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and 



H' 



1 



-3e 2 + 30 V3 1 ln2 \ )a 







e 2 \ln2\ 

as long as a is not too large. 

Of course the final judgement of the optimal c can be done by finding the minimum of H (c) with 
Mathematica or Matlab. The following are some examples. 



Graph of the modified MN curve withk=l(T 940 and 6=0.1 

kMN(c) 




131.06 131 



131.10 131.12 131.14 131.16 



Figure 6: Here n = 1,0 = — 1, a = 1, bo = 1 and S = 0.1. 
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Graph of the MN curve with 5=0.001 

MN(c) 




2 3 4 5 6* 

Figure 8: Here n = 1, (3 = -1, a = 1, b = 1 and 8 = 0.001. 
Graph of the MN curve with 5=0.0001 

MN(c) 




0.2 0.4 0.6 0.8 1.0 1.2 



Figure 9: Here n = 1,0 = -1, a = 1, & = 1 and S = 0.0001. 
Graph of the MN curve with 5=10" 5 

MN(c) 




0.2 0.4 0.6 0.8 1.0 1.2 

Figure 10: Here n = 1, /3 = -1, <7 = 1, b = 1 and 6 = 10~ 5 . 
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The often seen case f3 > is dealt with in the following Case3. 



Case3. j3 > and n > 1 Let / G E a and £ be the cube in Theoreml.3 with side length b . Let h 
be the map defined in (1) with /3 > and n > 1. For any fixed (5 in Theoreml.3, if < 5 < 47 (m+i) ' 
the optimal choice of c in the interval [12p v /ne 2 " 7 "7„(m + 1)5, oo) is the value c* minimizing -ff(c) 
on [12p v / ne 2 ™ 7 "7„(m + 1)5, oo), where 



(r) i±|±^ e cr 



#(c) := c 



<.ey 2 
e o 



ca + ^Jc 2 a 2 + 4<r(l + p + ■ 



Reason: This is a direct consequence of Corollary2.9. ft 

Remark: (i)If 1 + /3 — n > 0, -ff(c) is almost an increasing function of c. Therefore most time the 
smaller c is, the better it is in this situation. In other words, one should choose c = 12py / ne 2 ™ 7 ™7 Il (m+ 
1)5. This includes the most useful /? = 1 and n = 1 or 2. (ii)If 1 + /3 — n < 0, i?(c) -4- oo both as 
c — > + and c — > oo. However H(c) is diffcrcntiable and its minimum value does exist. In order to 
find the optimal c, one has to find where H'(c) — 0. 

Before going on to numerical results, let's analyze the critical point of H when 1 + (3 — n < 0. 
For 1 + P - n < 0, 



H(c) = 



1 + , l + <3 + n c{* -({*)-' 

c 4 (£ ) 4 e 2 e 2 " 



l+ff+T. 

= c 4 



ccr + ^cV 2 + 4<r(l + /3 + ra) 



> + Ve^ 2 + 4»(l+|i+~) 



T + Vc 2 g 2 +4g(l + /3 + n) 



G(c)e 



E(c) 



C 4 



where 



and 



G(c) := 



ca + ^Jc 2 a 2 + 4<t(1 + /3 + n) 



E(c) 

. Note that E(c) 

E'(c) = 



1 

2 

l 

32cr 



ccr + v /c 2 ct 2 + 4ct(1 + P + n) 



3ca - v / c 2 a 2 + 4cr(l + P + n) 
4^ 
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2c 2 cr 2 + 2co- v /c 2 o- 2 + 4cr(l + P + n) - 4cr(l + P + n) 

2c 2 a 2 

Aca + 2 v /c 2 cr 2 + 4cr(l + P + n) ' 



and 



^c 2 a 2 + 4a(l + /3 + n) 



> 



. Both G(c) and £7(c) are increasing functions of c. 
Now, 



ff'(c) 



[G'(c)e g ( c >+G(c)e g ( c )£'(c)] - G(c)e E ^c :l ^ 1 ^- 



C 2 



1G 



C 4 



cG'(c)e B ( c ) + C G(c)e E ( c )£'( C ) - n^l G (c)e E ^ 



n-l-fl 

C 2 



= 



iff 



. Note that 



n — 1 — 8 

cG'(c) + cG(c)E'(c) '-G{c) = 



da) = ^(d^ = ^^(n^-^cf 



and 



£ , (c) = I to-^lMH^n) 
V ; 2^ 4cr 



1 3ca-Vgg+4gI+g+g 1 

2 4cr c? 2 5 c 



3ccr - ^c 2 cr 2 + 4<r(l + /3 + n) 



where 



a 1 2a 2 c 

4 8 ~7?P + 4(7(1 + /3 + n) 



. Thus 
iff 



1 + / 3 + n c(r) ^4^-i 



ccr 2 + <7y/c 2 cr 2 + 4cr(l + /3 + n) 
4 v /c 2 cr 2 + 4cr(l + /3 + n) 

ff'(c) = 
ccr 2 + crVc 2 cr 2 +4cr(l + /3 + n) 



4 v /c 2 cr 2 + 4(7(1 + /? + 



n 



, . .... l+g+n 1 

+ c(£*) 4 2 

. . l±g±n 1 

+ c(D 4 2 
= 



3ccr - ^/c 2 <7 2 + 4cr(l + /3 + n) ca 2 + a^Jc 2 a 2 + 4cr(l + (3 + n) 



c 3 -,- 



4(7 4 v /c 2 cr 2 + 4ct(1 + + 

2ca 2 



4 8cr ^/cV 2 + 4a(l + /3 + n) 



iff 



l + /3 + n ,._i±4± i _ 1 



ccr 2 + cr(4£* - ccr) 



...... 1 

+ c(D 4 2 



n-l-/3 (r) i±^» 



4(4£* - ccr) 

3ccr - (4<f - ccr) ccr 2 + cr(4<f - ccr) 



4(7 



4(4£* - ccr) 



. 3 1 2cct 2 
+ C T- 



4 8cr (4£* - ccr) 



= 



17 



iff 



iff 



iff 



iff 



iff 
iff 
iff 



+ 



1 + p + n c 

4 e 



C<7 



16£* - 4ctr 



4ccr - 4<f ccr 2 + AaC 



4a 



16§* - 4ccr 



ccr 2 / 3 CCT 

+ ^ ^4 ~ 16£* - 



4cer 



n-l-p 



1+P + n c(4aC) c 
+ - 



4 

n - 1 - P 



(16f - 4ccr) 



(4c<7-4g*)(4<rg*) 
4cr 



+ T(12r -4ccr) 



= 



c(l + /3 + n)cr + ■ 



46cr(r) 2 + 16ctT 2 r 



+ 12(f) - 4cct£* 



4(7 

4(n - 1 - /3)£* + (n - 1 - /3)c<7 




c(l + p + n)aC + T ' 2 ' I-4(D 2 + 4ccrf + 12(D 2 - 4aj£*] 

- 4(n-l-/3)(D 2 + (n-l-/3)car 
= 



4c(f) 3 - 4(n - 1 - /3)(f) 2 + 2ncaC = 
2c(£* ) 2 - 2(n - 1 - /?)£* + nca = 

2(n - 1 - p) ± ^/4(n - 1 - pf - 



8nc 2 c 



4c 



(n - 1 - P) ± ^(n - 1 - /3) 2 - 2nc 2 er 



2c 



. Since there does exist c satisfying -ff'(c) = 0, it must be that 2nc 2 a < (n — 1 — /3) 2 , i.e. 



c < 



n-l-fi 



Numerical Result: Our numerical experience shows that c is usually near 



n-1- 



/2na 



Therefore 



n-l-(3 



is a good starting value when finding the optimal c iteratively by Mathematica. 

What's noteworthy is that H(c) increases rapidly if c > ■ Our experience shows that if c 

changes from to 10 (" 1 /3 ) or 10 °(" l-ff) ^. ne vame G r ( c ~) m ay expand 10 40 times. However 

-ff (c) expands very slowly as c goes from to 0+, although lim c ^ 0+ H(c) = oo. Of course, the 

requirement c e [12/9 v /ne 2 " 7 "7„(m + l)ci, oo) should be satisfied. Some examples are now offered. 
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MN(e) Graph of the MN curve with 6=0.0 1 




26.3 26.4 26.5 26.6 26.7 



Figure 11: Here n = 1,0 = 1, a = 1, b = 1 and 5 = 0.01. 
Graph of the MN curve with 5=0.001 

MN(c) 
500 '- 

400 - / 
300 1 / 
200 '- / 
100 - / 




Figure 13: Here n = 1, /3 = 1, a = 1, b = 1 and S = 0.0001. 
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Graph of the MN curve with 6=10 



MN(c) 



50 - 



40 - 



30 - 



20 - 



10 - 



Figure 14: Here n = 1, /3 = 1, a = 1, bo = 1 and S = 10 5 



Graph of the MN curve with 6=10 



MN(c) 



50 



40 



30 



2D 



10 



Figure 15: Here n = 1, /3 = 1, a = 1, 6 = 1 and 5 = lO" 6 . 



Now let's see an example which shows that n ^ n f ^ s a g°°d starting value for choosing the 



optimal c. In the following picture 



n-l-P _ 



0.408248. In fact a very simple Mathematica command 



shows that the critical point of the curve is exactly 0.408248. 
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Graph of the MN curve with 6=10 



MN(c) 




Figure 16: Here n = 3, /? = 1, a = 1, bo = 1 and (5 = 10 



-208 



4 How to choose c ?— a more theoretical approach 

In this section we take A^ of (4) into consideration. Theoretically < A < 1 and 5 can be arbitrarily 
small. Therefore A s is very influential. The value of A had been unknown for a long time. Fortunately 
it's clarified in [4]. This is a breakthrough and makes it possible to assess the influence of c on the 
error bound. However the value of A? highly depends on whether bo in Theoreml.3 is fixed or not. 
Hence we discuss it separately. 



4.1 b fixed 

Let bo in Theoreml.3 be fixed. Then 



C 



2p'^e 2n ^ ifce(0,co], 
3§- ifce[c ,oo) 



where Co := 3bopy/ne 2n7n . Since A = 

A* 



|) 6C " ln , we have 



e v(S)c if ce (0,co], 

(f) 1 ^ ifce[c ,oo) 



where V (S) := 



In* 



It's easily seen that A^ is a continuous function of c and is constant whenever c > c®. Its graph 
is as follows. 
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Graph of A* as function of c 



1 F 



0.00629300 k 




163.794 491.383 

Figure 17: Here n = l,/3 = -1,6 = M = 0.01 and < c < 491.382. 

A more interesting example can be seen in the following picture where c can be very large and 
5 is smaller than the preceding example. 

i 

Graph of A « as a function of c 




lxlO 204 2xl0 204 3xl0 204 4xl0 204 5 xl0 204 6xl0 204 



Figure 18: Here n = 3, /3 = -1, b = 1, 5 = 0.001 and 10" 5 < c < 6.15755 x 10 204 . 

Obviously is influential only when 8 is very small. The number j n grows very fast as n 
increases. Therefore A* ss 1 for high dimensions; unless 6 is extremely small. 
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By the graph of A^ one finds easily that if the optimal c lies on the right-hand side of Co when 
A 7 is not taken into consideration, it remains the same when A 7 is considered. 

With these understandings we can now begin our theoretical analysis of the optimal c. 



P = — 1 and n = 1 Let / € E a and E be the cube in Theoreml.3 with side length bo- 



Casel. 

Let h be the map defined in (1) with /3 = —1 and n = 1. For any fixed 6 in Theoreml.3, if 
< $ < 4 7 (m+i) 1 * ne optimal choice of c in the interval [12/9y / ne 2 ™ 7 "7„(rn + 1)6, 00) is as follows. 

(a) If Co < 12/9Y / ne 2 ™ 7 "7„(m + 1)5, choose c in the same way as Case2 of section3. 

(b) If Co > 12pY / ne 2n7 "7„(m + 1)6, the optimal c in the interval [12 p*Jne 2nin ^i n {m + 1)5, 00) 
is the number minimizing iJ(c)A? where was defined in Case2 of section3 and A^ is 
as in the beginning of subsection4.1. 

Reason: In (a) the optimal c lies to the right of cq. As explained in the beginning of subsectionLl, 
A* need not be considered and the choice of c is the same as section3. As for (b), the only difference 
between this section and section3 is that A^ should be considered. Hence we minimize H(c)\~x . jj 

All these can be done efficiently by Mathematica or Matlab. Note that the error bound grows to 
00 as c — y 00. The numerical process of finding the optimal c is about the same as Case2 of section3. 

Numerical Result:Let's see some examples. What's noteworthy is that for (b) the MN function 
is now iJ(c)A? . 

Graph of the modified MN Curve with k=10~ 945 and (5=0.1 

kMN(c) 



131.2 



131.4 



131.6 



131.8 



132.0 



Figure 19: Here n — 1,(3 — —1, a = 1, 60 = 1 and 5 = 0.1. 
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Graph of the MN Curve with 5=0.01 



MN(c) 



8x 10 



6x 10" 



4x 10 



2x10" - 




13.5 



14.0 



14.5 



15.0 



Figure 20: Here n = 1,P = -1, a = 1, b = 1 and S = 0.01. 
Graph of the MN Curve with (5=0.001 

MN(e) 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



Figure 21: Here n = 1,0 = -1, a = 1, b = l,and S = 0.001. 
^^Modified MN Curve with k=10 310 and 5=10 -4 

8xl0 302 



6xl0 302 - 



4x 10 302 



2xl0 302 - 




Figure 22: Here n = 1, /3 = -1, cr = 1, 6 = 1 and S = 10" 
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Graph of the modified MN Curve with k=10 830 and 6=IQ~ 5 

kMN(c) 



o.s 



0.6 



0.4 



0.2 



122 124 126 128 



130 



Figure 23: Here n — 1, (3 = —1, a = 1, &o = l,and 5 = 10 



Before going to the cases /3 > 0, n > 1 and ft = — 1, n > 2, some analytical work has to be done. 
Note that in both cases it satisfies 



ll/IU < fc c^ sup L II. 

[t&R™ e - J 

for any / g where fen is a constant independent of c. Moreover, 

Kl^eclCI _ | r |^e^*l 



sup 



iilff 

e " 



where | f | = ^yWp M, 

Therefore the crucial part of the error bound (4) is now 



H(c) 



I |f 



c ^ (iCl^e^-^Y ev(s)c if ce (0)Co]) 
c ^ (in^e^l-^) * (|) ^ if c e [co, co) 



(7) 



In order to use the results of section2, we assume that \n + f3\ > 1 and n + (3 + 1 > in the 
following case. 



Case2. 



\n + fi\ > 1 and n + (3 + 1 > 



Let f & Efj and -E be the cube in Theoreml.3 with side 
length bo. Let /i be the map defined in (1) with 1 + f3 + n > and \n + f3\ > 1. For any fixed 5 
in Theoreml.3, if < <5 < 47 (m+i) > * ne optimal c in [12 j o- v /ne 2 ™ 7 "7 rl (m + 1)<5, oo) is the number 
minimizing H(c) defined in (7). 
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Reason: This is just an immediate result of the paragraph preceding Case2. jj 

Remark: The most useful cases — — l,n > 2 and > 0,n > 1 are included here. In 
Case2, if 1 + — n > 0, then lim c _j. + ^( c ) = and lim c _ i . 00 ff(c) = oo. If 1 + — n < 0, 
lim c _ >0 + if (c) = lim c _>. 0o H(c) = oo. Finally, if 1 + — n = 0, lim c _ > g+ H(c) is a finite value and 
lim^oo H(c) — oo. Moreover, if 1 + ft — n = 0, H(c) is nearly increasing on [cq, oo). The optimal 
value of c can be obtained numerically by Mathematica or Matlab. 

Numerical Result:The following are some interesting examples. Note that in these examples 
only very smaMs are involved. Hence it's quite theoretical. We introduce the two most frequently 
seen cases n = 2, = — 1 and n = 1,0 = 1. 

Modified MN curve with k=10 69900 and 6=10~ 26 

kMN(c) 



8.xl(T 3 
6. x KT 3 

4.xitr 3 
2. x i(r 3 



1132 1134 1136 1138 1140 ° 

Figure 24: Here n = 2, = -1, a = 1, b = 1 and 6 = 1CT 26 . 
Graph of the modified MN curve with k=10 700 and <5=10" 25 

kMN(c) 
40 h 



30 



20 



10 



110 112 114 116 118 

Figure 25: Here n = 2, = -1, a = 1, b Q = 1 and 5 = 1CT 25 . 





2G 



Graph of the MN curve with 6=10 

MN(c) 




8 10 12 14 



Figure 26: Here n = 2,(3= — 1, a = 1, bo = 1 and <5 = 10 
Graph of the MN curve with <5=10~ 23 

MN(c) 



3 - 



2 - 



1 - 




2.0 2.5 3.0 3.5 4.0 4.5 5.0 




14.4 14.6 14.8 15.0 15.2 15.4 



Figure 28: Here n = 2, p = -1, a = 1, b = 1 and 8 = 10 
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MNM Graph of the MN curve with 5=0.0 1 



Ixl0 : 




26.25 26.30 26.35 26.40 26.45 



Figure 29: Here n = 1,0 = 1, a = 1, b = 1 and S = 10" 
Graph of the MN curve with 5=0.001 

MN(c) 




3.0 3.5 4.0 4.5 5.0 



Figure 30: Here n = l,/3 = l,a = I, bo = 1 and 8 = 10 
Graph of the modified MN curve with k=10 8 and (5=0.0001 

kMN(c) 




11.5 12.0 12.5 13.0 

Figure 31: Here n = 1,0 = 1,(7 = 1, bo = 1 and S = 10~ 
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Figure 32: Here n = l,/3 = l,a = 1, b = 1 and S = 1(T 5 . 
Graph of the modified MN curve with k=10 20550 and 5=W~ 6 

kMN(c) 




163.78 163.80 163.82 163.84 163.86 163.88 163.90 



Figure 33: Here n = 1,0 = l,a = 1, b = 1 and S = 1CT 6 . 



4.2 bo not fixed 

As explained in [5] and [8], some domains are invariant under dilation. Any point in such a 
domain is contained in a cube of side &o where bo can be made arbitrarily large and the cube is still 
contained in the domain. For example, 

fl := {(xi, • ■■,%): Xj > for i = 1, • • • , n} 

is such a domain. So is f2 = R n . 

In Theorem!. 3, if bo can be made arbitrarily large, then both C and A can be made arbitrarily 
small by increasing c and bo. The optimal choice of c will hence be very different. 

In this paper every approximated function / belongs to E a . The domain of / is of course R n . 
However the interpolation occurs in a cube as required in Theoreml.3. In practical problems the 
interpolation often can occur only in some subset f2 of R n , even if the domain of the approximated 
function is the entire R n . In this subsection the dilation- invariant domain f2 denotes the subset of 
R n where interpolation can occur. 

We begin with the case j3 = — 1 and n = 1. 
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Casel. P = — 1 and n = 1 Let f £ E a he the interpolated function and f2 C i? 1 be such that 
for any x € and &o > 0, there exists a cube i? of side bo such that a; G £ C O and interpolation 
can occur in E. Let /i be the map defined in (1) with /3 = — 1 and n — 1. For any 5 > 0, the optimal 
choice of c in the interval [12 / o A /ne 2 " 7 ™7 n (m + 1)5, 00) is the number minimizing H{c)e^ c where 
H{c) was defined in Case2 of section3 and rj(S) was defined in the beginning of sectionL 

Reason: The crucial part of the error bound in Theorem!. 3 is now H(c)\s where H(c) is de- 
fined as in Case2 of section3. Since C in Theoreml.3 can be kept equal to 2 / o'- v /ne 2 " 7 " by increasing 
bo, the value A^ will be equal to e^ 5 -^ for all c. The requirement c £ [I2p^/ne 2njn j n (m + 1)6, 00) is 
to make S < Sq in Theoreml.3. jj 



Remark: In Theoreml.3 &o appears first, and then Sq. However in the preceding case we reversed 
the order. We first fixed S, and then chose c. After c was chosen, we chose bg. The number b can be 
any real number greater than or equal to Zp ^ e 2 n ~, n ■ This will keep C = 2p'y / ne 2 ™ 7 " = 2 • ^ y/ne 2n " <n . 



Once bo was chosen, So 



} m+1 \ an d S < S automatically because c > 12/9 v /ne 2 ™ 7 "7„(m + 1)5, 



6C7„(rn+l) 



As a result, one should give S first and then choose c, and then 60. After bo is fixed, one can then 
arrange data points(interpolation points) in the cube E of side bo to make the fill distance d(E, X) 
equal to 5. Here d(E,X) := sup yeE mi x< £x \y — x\ as in Theoreml.3. 

What's noteworthy is that we only increase bo to keep C — 2//y / ne 2 " 7 ". We never decrease 60 
because it will only increase C and make A and So in Theoreml.3 worse. 



Numerical Result:Let's see some numerical examples as the following pictures show. The MN 
function is now H(c)e v ^ c . 



Graph of the modified MN curve with k=10" 1100 and 6=0.1 





kMN(c) 


l.xlO 


-166 




8.x 10 


-167 




6.x 10 


-167 




4.x 10 


-167 




2.x 10 


-167 





131.06 131.08 131.10 131.12 131.14 131.16 

Figure 34: Here n = 1, ^ = -1, a = 1 and S = 0.1 
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Graph of the MN curve with (5=0.01 

MN(c) 



1.5xl0 : 



1.0x10 



5.0x10 




13.2 13.4 13.6 13.8 14.0 14.2 14.4 



Figure 35: Here n = 1, = -1, a = 1 and 6 = 0.01 
Graph of the MN curve with 5=0.001 

MN(c) 




1.4 1.6 1.8 2.0 2.2 

Figure 36: Here n = 1,0 = -1, a = 1 and 6 = 0.001 
Graph of the MN curve with 5=0.0001 

MN(c) 




11 12 13 14 15 

Figure 37: Here n = 1,0 = -I, a = 1 and S = 0.0001 
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Graph of the modified MN curve with k=10 830 and (5=0.00001 

kMN(c) 



0.8 - 



0.6 



0.2 1 



122 



126 



I2S 



130 



Figure 38: Here n = 1,P = -1, a = 1 and 6 = 0.00001 



Let's now begin the second case. 



Case2. \n + @\ > 1 and n + /3 + 1 > Let h be defined as in (1) with |n + /3| > 1 and n + /3 + 1 > 0. 
Let / G E a be the interpolated function and il C R n be such that for any x G f2 and any 6o > 
there is a cube E of side 6o such that x G -E C Q and interpolation can occur in E. For any 6 > 0, 
the optimal choice of c in the interval [12pY / ne 2n7 "7„(m + 1)5, oo) is the number minimizing 



H(c) :-- 



\e 



where = ^ ^ 

was defined at the beginning of scction4. 

Reason: This case is different from Case2 of subsection4.1 only in A*. Here A^ is kept equal 
to e r '^ c by increasing bo- The criterion thus follows. 

Remark: As in Casel, 5 was given first. Then we chose c and then bo. The condition bo > 3p ^ e 2 n -, n 
is necessary. When applying this criterion, one should arrange data points to make X) = 5 
after 6q is chosen. 



Numerical Result: We begin our numerical examples from the case n = 2,0 = — 1 and then 
n = 1,13 = 1. 
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Graph of the modified MN curve with k=10 700 and <5=10" 25 

kMN(c) 



60 



40 



20 




Figure 39: Here n = 2, (3 = -1, a = 1 and 8 = 10 
Graph of the MN curve with <5=10~ 24 

MN(c) 



120 

-25 




Figure 40: Here n = 2, /3 = -1, a = 1 and 8 = lO" 24 . 
Graph of the MN curve with (5=10~ 23 

MN(c) 



30 : 

25 ; 
20 ; 
15 ; 
10 ; 

5 1 



Figure 41: Here n = 2,(3 = -l,cr = 1 and 5 = 10 23 . 
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Graph of the MN curve with 6=10 

MN(c) 




14.4 14.5 14.6 14.7 14.8 14.9 15.0 



Figure 42: Here n = 2, /3 = -1, a = 1 and S = 1CT 22 . 
Graph of the modified MN curve with k=l(T 1100 and c5=lCT 21 

kMN(c) 

2.0xl0 10 - / 
1.5 xlO 10 - / 
1.0x10'° - / 
5.0x10" - / 



142.92 142.94 142.96 142.98 143.00 143.02 



Figure 43: Here n = 2 )y 8 = -1, a = 1 and 6 = 1CT 21 . 
Graph of the modified MN curve with k=10" 4000 and 6=0.1 





kMN(c) 


1.2x10 


-269 




l.xlO 


-269 




8.x 10 


-270 




6.x 10 


-270 




4.x 10 


-270 




2.x 10~ 


-270 





262.075 262.080 262.085 262.090 262.095 

Figure 44: Here n = 1,0 = l,cr = 1 and 6 = 0.1. 
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MN(c) 
i x 10 40 - 



Graph of the MN curve with 5=0.01 



6xl0 40 



4x 10 41 



2x 10 40 



26.4 



26.6 



26.8 



27.0 



27.2 



Figure 45: Here n = 1, /3 = l,a = 1 and 6 = 0.01. 
Graph of the MN curve with 5=0.001 

MN(c) 
500 h 



4(10 



300 



200 



100 



Figure 46: Here n = 1, (3 = 1, a = 1 and 5 = 0.001. 
Graph of the MN curve with 5=0.0001 



MN(c-) 


1.x 10~ 7 




9.xl0- s 




8.X10- 8 




7.X10- 8 




6. x 10~ 8 




5. x 10- 8 




4.X10- 8 





12 



Figure 47: Here n = 1, 



13 14 15 

l,cr = 1 and 5 = 0.0001. 
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Graph of the modified MN curve with k=lCT u and (5=0.00001 

kMN(c) 




120 122 124 126 128 130 



Figure 48: Here n = 1, £ = 1, a = 1 and S = 0.00001. 
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